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Forecasting cosmological parameter constraints from near-future 

space-based galaxy surveys 

Anatoly Pavlov^, Lado Samushia^, and Bharat Ratra^ 
ABSTRACT 

The next generation of space-based galaxy surveys are expected to measure 
the growth rate of structure to about a percent level over a range of redshifts. 
The rate of growth of structure as a function of redshift depends on the behavior 
of dark energy and so can be used to constrain parameters of dark energy models. 
In this work we investigate how well these future data will be able to constrain 
the time dependence of the dark energy density. We consider parameterizations 
of the dark energy equation of state, such as XCDM and wCDM, as well as a 
consistent physical model of time-evolving scalar field dark energy, 0CDM. We 
show that if the standard, specially-flat cosmological model is taken as a fiducial 
model of the Universe, these near-future measurements of structure growth will 
be able to constrain the time-dependence of scalar field dark energy density to a 
precision of about 10%, which is almost an order of magnitude better than what 
can be achieved from a compilation of currently available data sets. 



1. Introduction 

Recent measurements of the apparent magnitude of Type la supernovae (SNela) con- 
tinue to i ndicate, quite convincingly, that the cosraological expans i on is currently acceleratiii g 
(see, e.g.. 



Conlev et al. II2OIII : ISuzuki et al. II2OI2I : iLi et al. II2OIII : iBarreira fc Avelinoll2011f ). 



If we assume that general relativity provides an adequate description of gravitational 
interactions on these cosmological length scales, then the kinematic properties of the Universe 
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can be derived by solving the Einstein equations 

1 



Here g^i, is the metric tensor, R^^, and R are the Ricci tensor and (curvature) scalar respec- 
tively, T^i, is the stress-energy tensor of the Universe's constituents, and G is the Newtonian 
gravitational constant. 

There is good observational evidence that the large-scale radiation and matter distri- 
butions are statistically spatially isotropic. The (Copernican) cosmological principle, which 
is also consistent with current observations, then indicates that the Friedmann-Lemaitre- 
Robertson- Walker (FLRW) models provide an adequate description of the spatially homo- 
geneous background cosmological model. 

In the FLRW models, the current accelerating cosmological expansion is a consequence 
of dark energy, the dominant, by far, term in the current cosmological energy budget. The 
dark energy density coul d be constant in time (and hence uniform in space) — Einstein's 



cosmological con stant A JPeebles 



varying in space (jPeebles &: Ratra 



9841) — or gradually decreasing in time and thus slowly 



19881). 



The "standard" model of cosmology is the spatially-flat ACDM model in which the cos- 
mological constant contributes around 75% of the current energy budget. Non-relativistic 
cold dark matter (CDM) is the next largest contributor, at around 20%, wit h non-relativistic 
baryo ns in third place with about 5%. For a review of the standard model see lRatra fc Vogeley 
(120081 ) and references therein. 



Recent m easurements of the anisotropies of the cosmic microwave background (CMB) 
radiation (e.g.. iKomatsu et al. 11201 ll : iReichardt et al. II2012I ). in conjunction with significant 
observationa l support for a low density of non-relativistic matter (CDM and baryons to- 



gether, e.g., IChen fc Ratral l2003l ) . as well as measurements of the p osition of the baryon 



Percival et al. 2010 



acoustic oscillation (BAO) peak in the matter power spectrum (e.g. 

Dantas et alTlboill : ICarnero et aTlboii lAnderson et al. Il2012h , provide significant observa- 
tional support to the spatially-fiat ACDM model. Other data are also not inconsistent with 
the standard ACDM model. These include strong grav itational lensing measurements (e.g., 
Chae et al. 112004 : iLee fc Ng]|2007l: iBiesiada et al. 1120101). measurement of Hubble parameter 



2011 



Brouzakis et al. 



2011 



as a function of redsh ift (e.g., Samushia fc Ratral 2006 : Sen fc Scherrei 20081: Pan et al. 



2010 



Chen fc Ratral 2011), large-scale structure data (e.g., Baldi fc Pettrino 2011 : De Boni et al. 



Carnpane 



Allen et al. 



2008 



li et al. 



20111). and galaxy cluster gas mass fr action 



measurements (e.g., 

recent reviews of the situation see, e.g., iBlanchardI (j2010l ). 
J201l[ ). 
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While the predictions of the ACDM model are in reasonable accord with current obser- 
vations, it is important to bear in mind that dark energy has not been directly detected (and 
neither has dark matter). Perhaps as a result of this, some feel that it is more reasonable to 
assume that the left hand side of Einstein's Eq. ([T]) needs to be modified (instead of postulat- 
ing a new, dark energy, component of the stress-energy tensor on the right hand side). While 
such modified gravity models are under active investigation, at present there is no compelling 
observational reason to prefer any of these over the standard ACDM cosmological model. 

The ACDM model assumes that dark energy is a cosmological constant with equation 
of state 

Pk = -pA, (2) 

where p\ and p\ are the pressure and energy density of the cosmological constant (fluid). This 
minimalistic model, despite being in good agreement with most observations available today, 
has some potential conceptual shortcomings that have prompted research into alternative 
explanations of the dark energy phenomenon]^ 

To describe possible time-dependence of the dark energy density, it has become popular 
to consider a more general equation of state parametrization 

= u{z)p^. (3) 

Here and are the pressure and energy density of the dark energy fluid with redshift 
z dependent equation of state parameter uj{z). The simplest such parametrization is the 
XCDM one for which the equation of state parameter is constant and results in accelerated 
expansion if oj{z) = ux < —1/3. In this case the dark energy density decreases with time 
and this allows for the possibility that the fundamental energy density scale for dark energy 
is set at high energy in the early Universe and the slow decrease of the energy density over 
the long age of the Universe ensures that the characteristic dark energy density scale now 
is small (a few meV). This also ensures that the dark energy density remains comparable 
to the matter energy density over a longer period of time (compared to that for the ACDM 
model). 

When LiJx = — 1 the XCDM parametrization reduces to the consistent (and complete) 
ACDM model. For any other value of ux the XCDM parametrization cannot consi stently 



descr i be spatial inhomogenei ties without further assumptions and extension (see, e.g.. iRatra 



1991uPodariu fc Ratral 120001 ). Models in which u{z) varies in time, cjCDM models, are also 



^ Structure formation in the ACDM model is governed by th e "standard" CDM structure formation 
mode l, which might be in some observational difficulty (see, e.g., Peebles fc Ratra 20031 : Perivolaropoulos 



20101). 
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unable to consistently describe spatial inhomogeneities without further assumptions and 
extension. 

A physically and observationally viable alternative to the ACDM model, that con- 
sistently describes a slowly decreasing in tirn e dark energy density, is the 0CDM model 
( jPeebles fc Ratralll988l : iRatra fc Peebles! 119881 ) . This model, in which a dark energy scalar 
field, slowly roles down its potential, resulting in a slowly decreasing dark energy density, 
alleviates some of the conceptual problems, mentioned above, associated with the ACDM 
model. The slowly rolling scalar field, at a given instant of time, can be approximated by a 
dark energy fiuid with an appropriately negative equation of state parameter. 

More specifically, a (pCDM model with an inverse-power-law scalar field potential energy 
density V^cj)) oc a > 0, is a prototypical example that has been extensively studied. 
This model has a non-linear attractor or "tracker" scalar field solution that forces the ini- 
tially sub-dominant dark energy density to come to dominate over the matter energy density, 
thus dominating the energy budget of the current Universe, and so resulting in the current 
accelerated cosmological expansion. In addition to therefore partially alleviating the "coin- 
cidence" problem of the ACDM model, the 0CDM model generates the current tiny dark 
energy scale of order an meV, measured by the SNela, through decrease, via cosmological 
expansion over the long age of the Universe, of a much larger energy scale. 

The a parameter controls the steepness of the scalar field potential, with larger val- 
ues resulting in a stronger time dependence of the approximate equation of state parameter 
and a = corresponds to t he ACDM model limit, a has been constrained using cur- 



rently available data (see e.g. , IChen fc Ratrall2004 IWilson et al. II2006I : IChen fc Ratrall2012 



Mania fc Ratrall2012l . and references t herein). The str ongest current limits are that a has 
to be less than ~ 0.7 at 2(T confidence (ISamushial l2009l ) . 



In the (f)CDM model, or in the XCDM or cuCDM parameterizations, the background 
evolution of the (spatially homogeneous) Universe differs from that in the ACDM case. This 
affects both the distance-redshift relation as well as the growth rate of large-scale structure. 
With precise measurements of distance and growth rate over a range of redshifts it will be 
possible to discriminate between cosmological models^! 

The BAO signature in the observed large-scale structure of the Univers e allows for the 
measurements of radial and angular distances as functions of redshift (see, e.g.. |Percival et al. 



^ There are many other models under current discussion, besides the ACDM and 0CDM models and 
XCDM and wCDM p arameterizat i ons w e consider here for illu s trative purposes. Fo r a sample of the avai l- 



able options see , e.g., Yang et al. I (l201ll). Frolov fc Guol (120111). iNunes et al. I (l201ll). iGrande et al. I (|2011[ ). 
Saitou fc NoUril |201ll) . Isilva et al. I |201ol ). iKamenshchik et al. I (|201l[ ). and lMaggiore et al. I |201ll) 
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2010l : iBlake et al. boilbl : ISeutler et al.boill : lAnderson et al. Il2012h . In addition, the redshift 
space distortion signal allows f or inferences about the strength of gravitational interactions 



on very large scale s fsee. e.g.. iPerciyal et al.l 120041: lAngulo et al. 1 120071 : iGuzzo et al.l 12008 



Blake et al. Il2011at ISamushia et al. II2012I : iReid et al. Il2012[ ). Currently available data sets 
have been used to measure distances and growth history up to a redshift z ~ 0.8 and the 
next generation of planned space-based galaxy redshift surveys of the whole extragalactic 
sky are expected to extend these measurements to a redshift 2; ~ 2. Possible candidates for 
such surveys i nclude the Euclid sat ellite mission that has been approved by the European 
Space Agency ( Laureijs et al. 201l[) and th e WFIRST satellite that was ranked high by the 



recent Decadal Survey ( IGreen et al. II2OIII ). These surveys have been shown to have the 
potential to measure angular distances, Hubble parameter H{z), and gr owth rate as func- 



tions of redshift to a few percent precision over a wide ra.nge o f redshifts (iWang et al.l 12010 
Samushia et al. II2OI1I : [Maierotto eralll2012HBasse et al.ir2012l )FI 



As mentioned above, an alternative potential explanation of the observed accelerated 
expansion of the Universe is to replace general relativity by a modified theory of gravity. For 
example, in the /(-R)-gravity models the Einstein-Hilbert gravitational action is modified to 



S 



16nG 



(4) 



where the function f{R) of the Ricci curvature R can in general be of any form. In the special 
case when f{R) = R one recovers the Einstein-Hilbert action which yields the Einstein equa- 
tions of general relativity, Eq. ([1]). For every dark energy model it is possible to find a function 



f(R) that will result in exactly the same expansion history (see, e.g. lSotiriou &: Faraomll2010 



Tsujikawall2010l : ICapazziello fc De Laurentisll201ll ) thus potentially eliminating the need for 
dark energy. However, nothing prevents the coexistence of modified gravity and dark energy, 
with both contributing to powering the current accelerated cosmological expansion. It is of 
significant importance to be able to determine which scenario best describes what is taking 
place in our Universe. 

In this paper we investigate how well anticipated data from the galaxy surveys mentioned 
above can constrain the time dependence of the dark energy. We will use the Fisher matrix 
formalism to obtain predictions for the 0CDM model and compare these with those made 
using the (model-dependent) XCDM and cuCDM parameterizations of dark energy. We will 
mostly assume that gravity is well described by general relativity, but will also look at some 
simple modified gravity cases. We find that the anticipated constraints on the parameter a 



^ For constraints on cosmological parameters from data from space missions proposed earlier, see 
Podariu et al. I (|200l[) and references therein. 
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of the (f)CDM model are almost an order of magnitude better than the ones that are currently 
available. 



Compared to the recent analysis of ISamushia et al. I (1201 ll ). here we use an updated 



characterization of planned next-generation space-based galaxy surveys, so our forecasts are 
a little more realistic. We also consider an add itional dark energy para metrization, XCDM, a 
special case of cuCDM that was considered by ISamushia et al. I (120111 ). as well as the (pCDM 
model, forecasting for which has not previously been done. 

The paper is organized as follows. In Sec. |2] we briefly describe the observables and 
their relationship to basic cosmological parameters. In Sec. |3] we describe the models of 
dark energy that we study. Section H] outlines the method we use for predicting parameter 
constraints, with some details given in the Appendix. We present our results in Sec. and 
conclude in Sec. O 



2. Measured power spectrum of galcLxies 

The large-scale structure of the Universe, which most likely originated as quantum- 
mechanical fluctuati ons of the scalar field that drove an early epoch of inflation (see, e.g.. 



Fischler et al. Ill985l ). became (electromagnetically) observable at z ~ 10^ after the recombi- 
nation epoch. Dark energy did not play a significant role at this early recombination epoch 
because of its low mass-density relative to the densities of ordinary and dark matter as well as 
that of radiation (neutrinos and photons). At 2 ~ 5 galaxy clusters began to form. Initially, 
in regions where the matter density was a bit higher than the average, space expanded a bit 
slower than average. Eventually the dark and ordinary matter reached a minimum density 
and the regions contracted. If an over-dense region was sufficiently large its baryonic matter 
collapsed into its dark-matter halo. The baryonic matter continued to contract even more 
due to its ability to lose thermal energy through the emission of electromagnetic radiation. 
This can not happen with dark matter since it does not emit significant electromagnetic 
radiation nor does it interact significantly (non-gravitationally) with baryonic matter. As 
a result the dark matter remained in the form of a spherical halo around the rest of the 
baryonic part of a galaxy. At z ~ 2 the rich clusters of galaxies were formed by gravity, 
which gathered near-by galaxies together. Also by this time the dark energy's energy density 
had become relatively large enough to affect the growth of large-scale structure. 

Different cosmological models with different sets of parameters can result in the same 
expansion history and so it impossible to distinguish between such models by using only 
expansion history measurements. This is one place where measurements of the growth 
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history of the large-scale structure of the Universe plays an important role. It is not possible 
to fix free parameters of two different cosmological models to give exactly the same expansion 
and growth histories simultaneously. It is therefore vital to observe both histories in order 
to obtain better constraints on parameters of a cosmological model. 

In a cosmological model described by the FLRW metric, and to lowest order in da rk mat 



ter o ver-density perturbations, the power-spectrum of observed galaxies is given by ([Kaiser 



19871) 



P,(A;,/i) = P„(A;)(6(T8 + M/i^)^ (5) 

Here subscript g denotes galaxies, Pm is the underlying matter power spectrum, h is the bias 
of galaxies, / is the growth rate, /i is the cosine of the angle between wave-vector k and 
the line-of-sight direction, and cxg is the overall normalization of the power spectrum (cxg is 
the rms energy density perturbation smoothed over spheres of radius 8/;,"^ Mpc, where h = 
i7o/(100kms^^Mpc~^) and Hq is the Hubble constant). Since, for a measured power spectrum 
of galaxies on a single redshift slice, the bias and growth rate are perfectly degenerate with 
the overall amplitude, in the equations below we will refer to bag and fag simply as b and /. 

The angular dependence of the power spectrum in Eq. can be used to infer the 
growth rate factor f{z) which is defined as the logarithmic derivative of the linear growth 
factor 

\ ^^^^ fR\ 

f{z) = — , (6) 

a ma 

where a is the cosmological scale factor, and the linear growth factor G{t) = 5(t)/5(tin) shows 
by how much the perturbations have grown since some initial time tinO 

The numerical value of the f{z) function depends both on the theory of gravity and on 
the expansion rate of the Universe. Since the growth rate depends very sensi tively on the 



total amount of non-relativistic matter, it is often parametrized as (see, e.g., iLinderl 12005 
and references therein) 

^ m^), (7) 

where 



and 



E^{z) ' 

E{z) = H{z)/Ho = + zy + + zy + noEiz). (9) 



* Here we have expanded the energy density p(t, x) in terms of a small spatially inhomogeneous fractional 
perturbation S{t,x) about a spatially-homogeneous background pb{t)'. p{t,x) ~ pb{t)[l + 5{t^x)\. 
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Here H{z) is the Hubble parameter and Hq is its value at the present epoch (the Hubble 
constant), Qrn is the value of the energy density parameter of non-relativistic matter at 
the present epoch {z = 0), Qk that of spatial curvature, and Qj:)e{z) is the energy density 
parameter which describes the evolution of the dark energy density and is different in different 
dark energy models. 

The growth index, 7, depends on both a model of dark energy as well as a theory of 
gravity. When general relativity is assumed and the equation of state of dark energy is taken 



to be of the general form in Eq. ([3]) then (see, e.g., lLinderll2005l . and references therein) 



7^ 0.55 + 0.05[l + a;(2 = 1)] (10) 

to a few percent accuracy. In the ACDM cosmological model 7 ~ 0.55. An observed 
significant deviation from this value of 7 will present a serious challenge for the standard 
cosmological model. 

The power spectrum is measured under the assumption of a fiducial cosmological model. 
If the angular and radial distances in the fiducial model differ from those in t he real cosmol- 



ogy, th e power spectrum will acqui re an additional angular dependence via the lAlcock fc Paczyhski 



(1 19791 . AP) effect, as discussed in ISamushia et al. I ( 1201 ll ) 



P^(,.,) . ( J- + - F-')) X |. + ^ _^^^ ]\ (11) 

where 

f\\{z)=Rr{z)/Rr{z), (12) 

f^iz) = DAiz)/DAiz), (13) 

F = fii/f±. (14) 

Here Rr = dr/dz is the derivative of the radial distance. Da is the angular diameter distance 
(both defined below), a hat indicates a quantity evaluated in the fiducial cosmological model, 
and a quantity without a hat is evaluated using the alternative cosmological model. The AP 
effect is an additional source of anisotropy in the measured power spectrum and allows for 
the derivation of stronger constraints on cosmological parameters. 
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3. Cosmological models 

In an FLRW model with only non-relativistic matter and dark energy the distances 
Da{z) and Rr{z) are 

1 / , dz' \ 

^^(^)= .,/n-n , ., sinhM4 / — - , (15) 



hy/rtk{l + Z) \ Jq E{z' 

Here E{z) is defined in Eq. (Q. The functional form of E{z) depends on the model of dark 
energy. 



3.1. ACDM, XCDM and wCDM parameterizations 

Here we describe the relevant features of the ACDM model and the dark energy param- 
eterizations we consider. 

If the dark energy is taken to be a fiuid its equation of state can be written as p = uj{z)p. 
For the ACDM model the equation of state parameter uj{z) = —1 and the dark energy density 
is time independent. 

In the XCDM parametrization uj{z) = ux{< ^1/3) is allowed to take any time- 
independent value, resulting in a time-dependent dark energy density. 

In the cuCDM parametrization the tim e dependence of coiz) is parametrized by intro- 
ducing an additional parameter Ua through (jChevallier fc Polarskil l200ll : iLinded |2003| ) 

W{z) = Wo + Wa——. (17) 
L -\- Z 

The XCDM parametrization is the limit of the wCDM parametrization with Ua = 0. In the 
ojCDM parametrization the function Qoe{z) that describes the time evolution of the dark 
energy density is 

and the corresponding expression for the XCDM case can be derived by setting Ua = here. 



^DEiz) = il-nm- + z)^^^+-o+«'.j exp I -3u;,^ ) , (1^ 
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3.2. (f)CBM model 

In the 0CDM model the energy density of the background, spatially homogeneous, 
scalar field (f) can be found by solving the set of simultaneous ordinary differential equations 
of motion, 

+ 3-0 + =0, (19) 
a 

8ttG , . k , . 

= ^(P + P0)--^, (20) 



1 ^1^2 + ^(^) ) , (21) 



Here an over-dot denotes a derivative with respect to time, a prime denotes one with respect 
to 0, V{(j)) is the potential energy density of the scalar field, is the energy density of the 
scalar field, and p that of the other constituents of the Universe. 



Following [Peebles fc Ratral ( 1l988h we consider a scalar field with inverse-power-law po- 



tential energy density 

V{<P) = ^0-". (22) 

Here a is a positive parameter of the model to be determined experimentally and k is a 
positive constant. This choice of potential has the interesting property that the scalar field 
solution is an attractor with an energy density that slowly comes to dominate over the 
energy density of the non-relativist ic matter (in the matter dominated epoch) and causes 
the cosmological expansion to accelerate. The function Qde^z) in the case of 0CDM is 

nDE{z) = ^{<P' + ^4>--). (23) 



4. Fisher matrix formalism 

The precision of the galaxy power spectrum measured in redshift bins depends on the 
cosmological model, the volume of the survey, and the distribution of galaxies within the ob- 
served volume. See App. A for a summary of how to estimate the precision of measurements 
from survey parameters. 

We assume that the power spectrum P{ki)'^^^^ has been measured in wave-number 
ki bins {i = 1 . . . N) and each measurement has a Gaussian uncertainty cTj. From these 
measurements a likelihood function 

C^exp(-^X'] (24) 
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can be constructed where 



X 



N 

E 

i=l 



Here p are the set of cosmological parameters on which the power spectrum depends. 



(25) 



The hkehhood function in Eq. (|24|) can be transformed into the hkehhood of theoretical 
parameters p by Taylor expanding it around the maximum and keeping terms of only second 
order in 5p as ^^{Sp) = Fj^SpPSp^^ where Fjk is the Fisher matrijo of the parameter set p 
given by second derivatives of the likelihood function through 

52 ln£' 



jk 



dpidp^ 



(26) 



The Fisher matrix predictions are exact in the limit where initial measurements as well 
as derived parameters are realizations of a Gaussian random variable. This would be the case 
if the P™^^ were perfectly Gaussian and the Pi{p) were linear functions of p, which would 
make the second order Taylor expansion of the likelihood around its best fit value exact. In 
reality, because of initial non-Gaussian contributions and nonlinear effects, the predictions of 
Fisher matrix analysis will be different (more optimistic) from what is achievable in practice. 
These differences are larger for strongly non-linear models and for the phase spaces in which 
the likelihood is non-negligible at some physical boundary (a = in case of 0CDM). A 
more realistic approach, that requires significantly more computational time and power, is 
to generate a large amount of mock data and perform a full Mqnte-C arlo Markov Chain 



(MCMC) analysis (see, e.g.. IPerotto et al. II2OO6 



Martinelli et al. 2011. where the authors 



find significant differences compared to the results of the Fisher matrix analysis). 

We assume that the full-sky space-based survey will observe Ha-emitter galaxies over 
1500 deg^ of the sky . For the density and b ias of observed galaxies we use predictions 



from lOrsi et al. I (I2OIOI ) and iGeach et al.l (120 lOl ) respectively. We further assume that about 
half of the galaxies will be detected with a reliable redshift. These numbers roughly mirror 
what proposed space missions, such as the ESA Euclid satellite and the NASA WFIRST 
mission, are anticipated to achieve. For the fiducial cosmology we use a spatially-flat ACDM 
model with VL^ = 0.25, the baryonic matter density parameter Qf, = 0.05, erg = 0.8, and the 
primordial density perturbation power spectral index Ug = 1.0, for convenience we summarize 
all the parameters of the fiducial model in Table 1. 



^ F or a review of the Fisher matrix formahsm as appUed to cosmological forecasting, see 
(|2n09h . 



Albrecht et al. 
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Table 1: Values of the parameters of the fiducial ACDM model and the survey. 









h 


0-8 


Us 


Efficiency 


Redshift span 


Covered sky area in deg^ 


0.25 


0.05 


0.0 


0.7 


0.8 


1.0 


0.45 


0.55 <z< 2.05 


15000 



We further assume that the shape of the power spectrum is known perfectly (for exam- 
ple from the results of the Planck satellite) and ignore derivatives of the real-space power 
spectrum with respect to cosmological parameters. 

We predict the precision of the measured galaxy power spectrum and then transform 
it into correlated error bars on the derived cosmological parameters. At first we make 
predictions for the basic quantities b and / in the XCDM and cuCDM parameterizations and 
in the 0CDM model. Then it allows us to predict constraints on deviations from general 
relativity and see how these results change with changing assumptions about dark energy. 
Finally, we forecast constraints on the basic cosmological parameters of dark energy models. 

For the XCDM parametrization these basic cosmological parameters are Pxcdm = 
{f,b, h,Qrn,^k,wx)- The cuCDM parametrization has one extra parameter describing the 
time evolution of the dark energy equation of state parameter, Pujcdm = (/, b, h, Qm, ^k, wq, Wa) 
For the 0CDM model the time dependence of the dark energy density depends only on one 
parameter a so we have p^cdm = {f,b,h,Qm,^k,C()- In order to derive constraints on the 
parameters of the considered cosmological models while altering assumptions about the cor- 
rectness of general relativity, we transform Fisher matrices of each model from the parameter 
set described above to the following parameter set (that now includes 7 that parametrizes the 
growth rate) Pmodei = (7, model), where by model we mean all the parameters of a particular 
model, for example, for cjCDM model = Pujcdm = (/, b, h, Q^, fi^, wq, Wa)- 



5. Results 

5.1. Constraints on growth rate 

Figure [1] shows predictions for the measurement of growth rate assuming different dark 
energy models. We find that in the most general case, when no assumption is made about 
the nature of dark energy, the growth rate can be constrained to a precision of better then 2% 
over a wide r ange of redshifts. This is in good agreement with previous similar studies (see. 



e.g.. Fig. 1 of ISamushia et al. II2OIII ). When we specify a dark energy model the constraints 
on growth rate improve by about a factor of two. There is very little difference between the 
results derived for different dark energy models: the precision is almost insensitive to the 
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assumed model. Also, one can notice that the curves for the XCDM parametrization and 
for the (pCDM model are almost identical. The likely explanation of this effect is that for 
a fixed redshift bin the 0CDM model is well described by the XCDM parametrization with 
the value of the parameter ux = P(t,/p<f>, where the values of the scalar field pressure and 
energy density are evaluated at that redshift bin. 

The measurements of growth rate can be remapped into constraints on parameters de- 
scribing the deviation from general relativity. Figure [2] shows correlated constraints between 
the current re-normalized Hubble constant h and the 7 parameter that describes the growth 
of structure. The 0CDM model constraints on both h and 7 are tighter than those for the 
XCDM or cuCDM parameterizations. As expected, the most restrictive ACDM model results 
in the tightest constraints. 



5.2. Constraints on dark energy model parameters 

We use measurements of growth and distance to constrain parameters of the dark energy 
models. 

Figure E] shows constraints on parameters of the cuC DM parametrization (these should 



be compared to Figs. 4a and 5a of ISamushia et al. II2OIII ) . When no assumptions are made 



about the nature of gravity the constraints on ojq and oja are very weak and degenerate. 
When we assume general relativity the constraints tighten significantly, resulting in ~ 10% 
accuracy in the measurement of uo and ~ 25% accuracy in the measurement of Ua- 

The upper panel of Fig. H] shows constraints on the parameters ux and Qm of the XCDM 
parametrization. Similar to the previous case, the constraints tighten significantly when we 
assume general relativity as the model of gravity. About a 2% measurement of ux and a 
5% measurement of Qm are possible in this case. The lower panel of Fig. Hlshow the related 
constraints on fi^ and Qm for the XCDM parametrization. The constraints are similar to, but 
somewhat tighter than, those for the wCDM parametrization. This is because the XCDM 
parametrization has one less parameter than the uCDM parametrization. Spatial curvature 
can be constrained to about 15% precision in this case. 

Figures 13 and show constraints on parameters of the 0CDM model. In the most general 
case, when no assumption is made about the nature of gravity, the constraints are weak and 
the parameters a and Qm are strongly correlated, with larger values of a requiring larger 
values of Qm- When general relativity is assumed, the constraints become much stronger 
and parameter a can be constrained to be less than 0.1 at the 1-a confidence level. This is 
significantly better than any constraint available at the moment. 
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Figure [7] shows constraints on the parameters of the ACDM modeL From the clustering 
data alone the spatial curvature can be constrained with almost 1% precision, largely because 
this model has the least number of free parameters. 

The exact numerical values for the forecast error bars and likelihood contours should be 
taken with caution and not be interpreted as predictions for the performance of any specific 
survey (such as Euclid or WFIRST). Our main objective in this work was first to investigate 
how the modified gravity constraints change with different models of dark energy and second 
to demonstrate the improvement in 0CDM model constraints achievable with future galaxy 
surveys. Because of this we were able to simplify our method by adopting a Fisher matrix 
formalism instead of a full MCMC approach and also use a simplified des cription of the 



survey baseline. For mor e realis t ic predictions of E ' ttc/id p erformance, see, e.g.. lLaureijs et al. 
feoilh : [Samushia et al. I fl201ll ): iMaierotto et al.1 f |2012h . 



6. Conclusion 

We have forecast the precision at which planned near-future space-based spectroscopic 
galaxy surveys should be able to constrain the time dependence of dark energy density. For 
the first time, we have used a consistent physical model of time-evolving dark energy, (/)CDM, 
in which a minimally-coupled scalar field slowly rolls down its self-interaction potential energy 
density. We have shown that if general relativity is assumed, the deviation of the parameter 
a of the 0CDM model can be constrained to better than 0.05; this is almost an order of 
magnitude better than the best currently available result. 

The constraints on basic cosmological parameters, such as the relative energy densities 
of non-relativistic matter and spatial curvature, depend on the adopted dark energy model. 
We have shown that in the 0CDM model the expected constraints are more restrictive than 
those derived using the XCDM or wCDM parameterizations. This is due to the fact that 
the 0CDM model has fewer parameters. Also, the XCDM and wCDM parameterizations 
assign equal weight to all possible values of cu, while in the </)CDM model there is an implicit 
theoretical prior on which equation of state parameter values are more likely, based on how 
easy it is to produce such a value within the model. 

Since the observational consequences of dark energy and modified gravity are partially 
degenerate, constraints on modified gravity parameters will depend on the assumptions made 
about dark energy. In Table [3] we show how the constraint on the 7 parameter depends on 
the adopted dark energy model. The constraints on 7 are most restrictive in the ACDM 
model. For the 0CDM model the constraints on 7 are about a third tighter than those for 
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the cuCDM and XCDM parameterizations. 

These results are very encouraging: data from an experiment of the type we have 
modeled will be able to provide very good, and probably revolutionary, constraints on the 
time evolution of dark energy. 
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A. Appendix 

In this Appendix we summarize how to estimate the precision of measurements from 
the survey parameters. 

The Fisher matrix coefficients are given by 
where the effective volume is 

_ nP{k,fi) 

- ^\ + nP{k,f,y ^^^^ 



and Vq is the total survey volume and n is the number density. Also, following iTegmark 



(Il997l ). we multiply the integrand in Eq. flAip by a Gaussian factor exp ( — fc^a^^^^^ ) , where 



dz 

r{z) is the co- moving distance, in order to account for the errors in distance induced by the 
errors of redshift measuremen ts, = 0.001. We model the theoretical power spectrum using 



an analytic approximation of lEisenstein &: Hul ( 119981 ). We integrate in k from fcmin = to 



^max5 where the /Cmax values depend on redshift and are chosen in such a way that the small 
scales that are dominated by non-linear effects are excluded. The range of scales that will 
be fitted to the future surveys will depend on how well the theoretical templates are able 
to describe small-scale clustering and is difficult to predict. The A;max values along with the 
expected bias and number density of galaxies are listed in Table [2l 



In order to derive the Fisher matrix of a specific cosmological model we have to go from 
our initial parameter space to the parameter space of the cosmological model whose Fisher 
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matrix we want. Th e transformation formula for the Fisher matrix is given by (see, e.g., 



Albrecht et ah 2009. for a review) 



where the primes denote the "new" Fisher matrix and parameters. 

We now hst the derivatives of the transformation coefficients of the 0CDM model in the 
limit a — > and flk — > (which corresponds to the fiducial spatially-flat ACDM model). 
The transformation coefficients relating f\\{z) and the parameters {h,Qrn,^k,Ci) are 



df\\{z) _ 1 
dh h 

dfiiiz) _ 1 

dQrr, 2E(z 



(A4) 

-,[l-il + z)% (A5) 



= - (I + z?] 

df\\{z) _ {1-nj 



da 8E{zy 

For the other transformation coefficients, it is convenient to introduce the integral 



(A6) 
(A7) 



^^'^-Iw-y '''' 

Then the transformation coefficients between f±{z) and the parameters {h,Qm,^k, «) are 

i-{i+zn (Aio) 



dhiz) 1 r dz' ,^ , 



d^^ 2D{z)Jo E{z'f 
df^{z) D{zf 1 r dz' 



+ / _;i^[i-(i + /)2i (All) 

(A12) 



90fe 6 2D{z) Jq E{z' 

df^iz) (l-Qm) r dz' 



da 8D{z) Jq E{z'f' 

Finally, the transformation coefficients between the growth factor f{z) and the param- 
eters (7, /i, f^fe, a) are 

dfjz) f{z) 

-^ = —-\nf{z), (A13) 
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dfjz) _ ^f{z) 



{Eizy ~ + zf - 1]} 



^ = -^\(l + z?-l] 

on, E(z)^^^^ ' J' 



dfjz) 
da 



iEizV 



(A14) 
(A15) 
(A16) 



Table 2: Va lues of the /cmn-r, bia s b{z) 
taken from lGeach et al.l (j201Cll ). 



from lOrsi et al. I ( l2010l ). and the number densities n{z) 



2, 


k 

"'max 






0.55 


0.144 


1.0423 


3220 


0.65 


0.153 


1.0668 


3821 


0.75 


0.163 


1.1084 


4364 


0.85 


0.174 


1.1145 


4835 


0.95 


0.185 


1.1107 


5255 


1.05 


0.197 


1.1652 


5631 


1.15 


0.2 


1.2262 


5972 


1.25 


0.2 


1.2769 


6290 


1.35 


0.2 


1.2960 


6054 


1.45 


0.2 


1.3159 


4985 


1.55 


0.2 


1.4416 


4119 


1.65 


0.2 


1.4915 


3343 


1.75 


0.2 


1.4873 


2666 


1.85 


0.2 


1.5332 


2090 


1.95 


0.2 


1.5705 


1613 


2.05 


0.2 


1.6277 


1224 
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Fig. 1. — Predicted relative error on the measurements of growth rate as a function of 
redshift z in redshift bins of = 0.1 for different models of dark energy. The upper solid 
black line shows predictions for the case when no assumption is made about the nature of 
dark energy. 
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Fig. 2. — Predicted one standard deviation confidence level contour constraints on the current 
renormalized Hubble constant h and the parameter 7 that describes deviations from general 
relativity for different dark energy models. 
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Fig. 3. — Upper panel shows one standard deviation confidence level contours constraints 
on parameters oUa and uuq of the a;CDM parametrization, while lower panel shows these for 
parameters and Q^. 
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Fig. 4. — One standard deviation confidence level contour constraints on parameters of the 
XCDM parametrization. 
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Fig. 5. — One standard deviation confidence level contour constraints on parameters a and 
of the 0CDM model. Lower panel shows a magnification of the tightest two contours in 
the lower left corner of the upper panel. 
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Fig. 6. — One standard deviation confidence level contour constraints on parameters and 
Qjn of the 0CDM model. The lower panel shows a magnification of the two tightest contours 
in the center of the upper panel. 
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Fig. 7. — One standard deviation confidence level contour constraints on parameters and 
flm of the ACDM model. 



Table 3: Predicted deviations of parameter 7 from its fiducial value, at one standard deviation 
confidence level, for different assumptions about dark energy. 



DE model 


Fiducial 7 


deviation 


a;CDM 


0.55 


0.035 


0CDM 


0.55 


0.023 


XCDM 


0.55 


0.035 


ACDM 


0.55 


0.016 




